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ABSTRACT 

We constrain the shape of the primordial power spectrum using recent measurements of 
the cosmic microwave background (CMB) from the Wilkinson Microwave Anisotropy 
Probe (WMAP) 7-year data and other high-resolution CMB experiments. We also in- 
clude observations of the matter power spectrum from the luminous red galaxy (LRG) 
subset DR7 of the Sloan Digital Sky Survey (SDSS). We consider two different mod- 
els of the primordial power spectrum. The first is the standard nearly scale-invariant 
spectrum in the form of a generalised power-law parameterised in terms of the spectral 
amplitude A s , the spectral index n s and (possibly) the running parameter n run . The 
second spectrum is derived from the Lasenby and Doran (LD) model. The LD model 
is based on the restriction of the total conformal time available in a closed Universe 
and the predicted primordial power spectrum depends upon just two parameters. An 
important feature of the LD spectrum is that it naturally incorporates an exponential 
fall-off on large scales, which might provide a possible explanation for the lower-than- 
expected power observed at low multipoles in the CMB. In addition to parameter 
estimation, we compare both models using Bayesian model selection. We find there 
is a significant preference for the LD model over a simple power-law spectrum for a 
CMB-only dataset, and over models with an equal number of parameters for all the 
datasets considered. 
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- cosmic background radiation - large-scale structure 



1 INTRODUCTION 

Cosmological inflation not only explains the homogeneity 
of the universe on large scales, but also provides a theory 
for explaining the observed level of anisotropy ( Guth||1981 



Albrecht fc Steinhardt||1982| |Mukhanov fc Chibisov||i982 i 



Inflationary models generically predict Gaussian, adiabatic 
and nearly scale-invariant primordial fluctuations. To deter- 
mine the shape of the primordial power spectrum from cos- 
mological observations, it is usual to assume a parameterised 
form for it. The simplest assumption is that the initial spec- 
trum has the form of a simple power-law, parameterised in 
terms of the spectral amplitude A s and the spectral index 
n s . Recent analyses have shown, however, that the spectral 
index may deviate from a constant value (close to unity) and 
so consideration of models that provide some running of the 
index (defined by n run = dn s /d\nk) is warranted. Further- 
more, models that predict a decrement in CMB power at low 



multipoles seem to be preferred by observations (Efstathiou 



2003 Contaldi et al.|2003 l, so that, an extra variable to ac- 
count for a cut-off scale might be considered. It would be 
feasible to continue adding parameters in this fashion un- 
til some arbitrary accuracy of model fit is achieved, but this 
procedure would fail to account for Occam's razor: a simpler 
model should be preferred, unless the data require a more 
sophisticated one. Using the Bayesian evidence to select be- 
tween models is one way to include this consideration. 

There have been several recent studies regarding the 
shape of the primordial spectrum, some based on physical 
models, some using observational data to constrain an a pri- 
ori parameterisation, and others attempting a direct recon- 
struction from data ( Barriga et al.|2001| |Bridle et al.|2003| 
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A primary goal of this work is to fit the predicted 
form of the primordial spectrum from a physically moti- 
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vated model described by Lasenby & Doran ( 2005 1 (hence- 
forth LD). We translate observational data into constraints 
on the generalised power-law and LD forms for the primor- 
dial power spectrum (and the standard cosmological param- 
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eters), and decide which model provides the best fit to ob- 
servational data using the Bayesian evidence. 

The paper is organised as follows: in Section 2 we 
present the two different models for the power spectrum 
and in Section 3 we describe basic parameter estimation 
and model selection. We list the datasets and the cosmo- 
logical parameters considered in Section 4 and present the 
resulting parameter constraints in Section 5. We compute 
Bayesian evidences in Section 6 to decide which model pro- 
vides the best description for current observational data and 
we validate our analysis by applying it to a CMB simulated 
data. Our conclusions are presented in Section 7. 



2 PRIMORDIAL POWER SPECTRUM 

The correlation function £ of density fluctuations S = Sp/p 
at two separated points x and x + r is defined as 



£(r) = (S(x)5(x + r)>. 



(1) 



Because the assumption of homogeneity and isotropy, £ is a 
function only of r = jrj. The power spectrum P(k) describes 
the amplitude of fluctuations on different length scales and 
it is related with the inverse Fourier transform of the corre- 
lation function £ by: 



P(k) 



(|4| 2 ). 



(2) 



During the inflationary period, fluctuations in the inflaton 
field S(j>k result in curvature perturbations IZ(k) given by 



TZ(k) = - 



H 



(3) 



where the quantities are evaluated at the horizon exit epoch 
k — RH. Here R is the scale factor of the universe and 
H = R/R is the Hubble parameter. In this paper, we follow 
a slow-roll approximation (e.g. Liddle fc Lyth|1999 |. Hence 
the power spectrum of the inflaton fluctuations is constant 
in time and equal to 



Ps4,(k) = 



(4) 



Thus, the primordial curvature spectrum Piz(k) computed 
from |2| - Q is 



Vn(k) = 



H 



(5) 



index to vary as a function of scale, such that n s (k). This 
can be achieved by including a second order term in the 
expansion of the power spectrum 

/ h \ ,I s _1 + ( 1 / 2 ) ln(k / kg) (dn/d In k) 

where dn/dlnk is termed the running parameter n Iun and 
we would expect n run as for standard inflationary models. 

In what follows we will consider three Power-Law pa- 
rameterisations. In the first model (PL 1), we will assume a 
simple power-law spectrum (without running) and restrict 
the universe to be spatially flat. In PL model 2, we allow the 
spatial curvature of the universe to be a free parameter. In 
PL model 3, we allow for a running spectrum, but again re- 
strict the universe to be spatially-flat. In this way, power-law 
models 2 and 3 have the same number of parameters. 

2.2 The LD model 

Assuming the cosmological constant is the origin of dark en- 



ergy, Lasenby & Doran ( 2005 I provided a construction for 



embedding closed- universe models in a de Sitter background. 
As a consequence of this novel approach, a boundary condi- 
tion on the total available conformal time emerges. Defining 
the total conformal time r\ as 

f°° dt 

the LD model requires r\ — it/2. For more details about the 
choice of the boundary condition, including how it can be 
reinterpreted as an eigenvalue condition on the solution of a 



differential equation, see Lasenby ( 2003 1 ; Lasenby & Doran 
( |2004| |2005[ ) . In order to understand some consequences of 
the new boundary condition we split the history of the Uni- 
verse in two main contributions to the total conformal time: 
matter (radiation and dust) and inflationary eras. Hence, we 
want to compute the conformal time r/M elapsed during the 
matter era and add it to that elapsed in the inflationary era 
rji, such that the boundary condition is satisfied: 



7)i + Vm = 



(9) 



It is found that this constraint leads to a 'see-saw' mecha- 
nism linking the parameters describing the current state of 



the universe with the initial conditions (Lasenby & Doran 
[20041). 



2.1 Power-law parameterisations 

Cosmological slow-roll inflation predicts the spectrum of 
curvature perturbations to be close to scale-invariant. Based 
on this, the spectrum is commonly assumed to have the form 

/ b \ ns ~ 1 

Ms) - (6) 

where the spectral index n s is expected to be close to unity; 
ko is the pivot scale (set to ko = 0.05 Mpc -1 throughout). 
A spectrum where the typical amplitude of perturbations is 
identical on all length scales is known as Harrison-Zel'dovich 
spectrum (n s = l).This particular parameterisation involves 
only one free parameter, the spectral amplitude V(k) = A s . 
A further extension is possible by allowing the spectral 



2.2.1 Matter era 

The general description of the large scale Universe is based 
on the Robertson- Walker space-time with dynamics gov- 
erned by the Einstein equations. The resulting Friedmann 
equations can be written as (with c = 1) 





R 2 + k 
R 2 


A 
3" 


8tvG 


4 


R 2 + k 
+ R 2 


A = 


SttGP. 



(10) 



Here k — 0, ±1 defines the geometry of the universe, A is the 
cosmological constant, and the relationship between density 
p and pressure P is encoded in the equation of state P = 
7p. The behaviour of the homogeneous universe is governed 
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by the parameters representing its matter-energy content, 
namely 



ft M = 



8nGp 
3H 2 ' 



A = 



A 



(ii) 



and its expansion history defined by the Hubble parameter 
H. Moreover, if we assume the matter density is made up of 
decoupled dust and radiation, the equations governing Qm 
and Qa can be solved exactly and its solution is controlled 
by two arbitrary constants a and /3 given by 



(fir, 



a/3 



1)3 : 



(12) 



(n mo + n ro + n Ao - i) 2 



where subscript 'o' denotes quantities evaluated at present 
time. The total conformal time for this type of universe can 
be written in terms of the dimensionless parameters a and 
f3 as (see |Lasenby &: Doran|2005| ) 



Vm = 



f 

Jo 



dx 

(fix 4 + x 3 -x 2 + a y/ 2 ' 



(13) 



2.2.2 Inflationary era 

The computation of the conformal time in the inflationary 
epoch is a more elaborate process. Let us consider a basic 
inflationary model where the particle responsible for this 
process is simply a real, time-dependent, homogeneous, free, 
massive scalar field cj>, described by the equations 



H + H--- + 



A 4vrG 



(20 2 - mV) = 0, 



■ + 3H(j> + m 2 



0. 



(14) 



For closed universe models, the scale factor is given explicitly 
by 



1 47rG.- 2 2 2 

W = — (0 +m0) 



„ 2 A 
~3' 



(15) 



In order to compute the conformal time rji , it is necessary to 
seek out suitable conditions before the onset of inflation and 
then solve the dynamics for the scalar field encoded in equa- 



tions (14 1. To do this, Lasenby & Doran (20051 developed 
formal series expansions out of the initial singularity, t = 0, 
in terms of dimensionless variables u — t/t p and fj, — m/m p , 
where the subscript 'p' denotes a Planck units variable. The 
series are given by 



4>{u) 



1 00 



H(u) = ijTff„(u)ln>), 



(16) 



with 



b () + 6 4 u 4/3 
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(17) 



3 3u 3 

We observe that two new free parameters 6o and 64, 
appear in the series expansions in \17\ . Together with the 
mass of the scalar field /i, they control the magnitude of 
the field and how long the inflationary period lasts. In order 
to decide on the priors we shall employ in our subsequent 
Bayesian analysis, it is worth pointing out some features 
related with these new parameters. 

• The amplitude of the perturbations is determined by 
the scalar field mass [i. To match the observed level of CMB 
anisotropics, we shall need to set it to be about 

-6 



l«~ 10' 



(18) 



• The number of e-foldings N is primarily determined by 
bo and may be approximated as 



N : 



27T&0. 



(19) 



Hence, to obtain realistic models we need bo to be of order 
~ a few. 

• The conformal time is estimated by 



Employing the constraint |?>4|/i~ 4 / 3 should thus be 

around unity. 

• The parameter 64 controls the initial curvature, as can 
be seen from (Tl5l): 



R 

L 



( 2187 \ 
V 12544tt J 



V* „l/3 



(21) 



Therefore 64 must be negative. Making use of the rest of the 
parameters and (201, |&4| should be around 10 -9 . 



The restriction on the values for the model parameters to- 
gether with the boundary condition, severely limits the class 
of models allowed to reproduce current cosmological obser- 
vations. 

Once we have found the initial conditions (171, it is 



straightforward to solve numerically the dynamics of the 
scalar field <f> and the expansion history H to determine 
the evolution of the universe. As an example, let us con- 
sider the best-fit values for the cosmological parameters 



given by WMAP+BAO+#o (Komatsu et al. 2011 1, in which 



case the conformal time elapsed during the matter epoch is 
T]M = 0.22. In order to satisfy the boundary condition we 
should choose appropriate values {bo, 64, /i} such that by the 
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Table 1. Jeffreys scale for evaluating the strength of evidence 
when two models are compared. 

| In B i| Odds Probability Strength of Z 

< 1.0 < 3 : 1 < 0.750 Inconclusive 

1.0 ~ 3 : 1 0.750 Significant 

2.5 ~ 12 : 1 0.923 Strong 

> 5.0 > 150 : 1 > 0.993 Decisive 



ln(u) 



Figure 1. Evolution of the conformal time ni as a function of 
ln(u). We observe that r\i saturates at a value of around 1.35 
by the end of inflation. The parameters used in this model are 
b = 2.47, f> 4 = -17.7 X 10~ 9 and u = 1.68 X 10" 6 . 



end of the inflationary period the achieve conformal time is 
rji ~ 1.35, as shown on Fig [T] The expected shape of the 
primordial power spectrum V(k) for this model is directly 
computed from ([sj . 



3 BAYESIAN ANALYSIS 
3.1 Parameter Estimation 

A Bayesian analysis provides a consistent approach to es- 
timate the set of parameters © within a model M, which 
best describes the data D. The method is based on the as- 
signment of probabilities to all quantities of interest, and 
then the manipulation of those probabilities given a series 
of rules, of which Bayes' theorem plays the principal role, 
(Baycs 1761). Bayes' theorem states that 

P(D\&,M) P(&\M) 
P{D\M) ~~ ' 

where the prior probability P(&\M) = tt represents our 
knowledge of before considering the data. This probabil- 
ity is modified through the likelihood P(D\&, M) = C when 
the experimental data D are considered. The final objective 
for Bayesian inference is to obtain the posterior probability 
P(&\D, M) which represents the state of knowledge once 
we have taken the new information into account. The nor- 
malisation constant is the marginal likelihood or Bayesian 
evidence P(D\M) = Z. Since this quantity is independent 
of the parameters 0, it is commonly ignored in parame- 
ter estimation but takes the central role for model selection 
( |Hobson et al.|2002||Llddle|2004) . 



P(@\D,M) = 



(22) 



3.2 Model selection 

A complex model that explains the data slightly better than 
a simple one should be penalised for the introduction of extra 
parameters, because the additional parameters bring with 
them a lack of predictability. Also, if a model is too simple, it 
might not fit the data equally well, then it can be discarded 



(e.g. |Liddle|[2006| |2007| |Trotta|[2008[ ). Following this line, 
many attempts have been performed to translate Occam's 
razor into a mathematical language for model selection. The 



prime tool for the model selection we focus on is the Bayesian 
evidence. 

The Bayesian evidence provides a natural mechanism 
to balance the complexity of cosmological models and then, 
elegantly, incorporates Occam's razor. It applies the same 
type of analysis from parameter estimation but now at the 



level of models rather than parameters (Trotta 20071. Let 



us consider several models M, each of them with probability 
P(M). The Bayesian evidence appears again, but now in 
another form of the Bayes theorem: 



P(M\D) 



P(D\M)P(M) 
—pjD) 



(23) 



The left-hand side denotes the probability of the model given 
the data, which is exactly what we are looking for model 
selection. 

It was previously mentioned that the Bayesian evidence 
is simply the normalisation over its posterior expressed by: 



Z = J C(@)n(@)d N Q, 



(24) 



where N is the dimensionality of the parameter space. When 
two models are compared, Mo and Mi , the quantity to bear 
in mind is the ratio of the posterior probabilities given by 



P(M \D) _ Z P(Mo) 



P{Mx\D) 2iP{Mx) 



Br 



P(M ) 
P(Mi) 



(25) 



where P{Mq)/P(M\) is the prior probability ratio for two 
models, usually set to unity. The evidence ratio, often 
termed as the Bayes factor Boi , quantifies how well model 
may fit data when it is compared to model 1. Jeffreys| ( |1961[ ) 
provided a useful guide on which we are able to make qual- 
itative conclusions based on this difference (see Table [T]). 

Until recently, numerical methods such as thermody- 
namic integration ( Beltran et al.| |2005 Bri dges et al.|[2007[ ) 
required around 10' likelihood evaluations to obtain accu- 
rate estimates of the Bayesian evidence, making the pro- 
cedure hardly tractable. The nested sampling algorithm, in- 



vented by Skilling ( 2004 ) , has been successfully implemented 



for cosmological applications ( Mukherjee et al.]|2006| |Shaw| 
et al. 1120071 iFeroz et al. 2008 I requiring about a hundred 



times fewer posterior evaluations than thermodynamic in- 
tegration to achieve the same accuracy in the evidence es- 
timate. A substantially improved and fully-parallelized al- 
gorithm called MultiNest ( |Feroz et al.|2009 1 increases the 
sampling efficiency for calculating the evidence and obtain- 
ing posterior samples, even from distributions with multi- 
ple modes and/or pronounced degeneracies between param- 
eters. In our case, to carry out the exploration of the cosmo- 
logical parameter space we use a modified version of both 
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the MultiNes t ( |Feroz et a1"1|2009| ) and CosmoMC ( |Lewis 
fc Bridle||2002[ ) packages. 



Table 2. Parameters and prior ranges used in our analysis. 



4 DATASETS 

Measurements of the CMB anisotropies and large-scale 
structure play an important role in both fitting parame- 
ters and comparison of models in cosmology. To constrain 
the space-parameter in each model, we first use the latest 



7-year data release from WMAP (henceforth WMAP7; Lar- 
son et ah|2011 1, which includes a good measurement up to 
the third acoustic peak in the temperature CMB spectrum. 
This comprises our dataset 1. 

WMAP7 data by itself cannot, however, place strong 
constraints on all the parameters because of the existence 
of parameter degeneracies, such as the r — A degeneracy 
and the well-known geometrical degeneracy, involving f2 m , 
Qa and Qk- Nevertheless, when WMAP7 is combined with 
other cosmological observations, they together increase the 
constraining power and considerably weaken the degenera- 
cies. In addition to WMAP7, we therefore include recent 
results from CMB experiments that are able to reach higher 
resolution on small patches of the sky, such as the Arcminute 
Cosmology Bolometer Array (AC BAR; |Kuo et al 



2004). 



Cosmic Background Imager (CBI; Readhead et al.||2004[ ), 
Ballon Observations of Millimetric Extra-galactic Radiation 
and Geophysics (BOOMERang; | Jones et al.|2006| ). We also 
include observations of the matter power spectrum from the 
luminous red galaxy (LRG) subset DR7 of the Sloan Digital 
Sky Survey (SDSS; |Reid et al.|[20To| . In addition to CMB 
and galaxy surveys we include the Hubble Space Telescope 



(HST; |Riess et aL]|2009[ ) key project for the Hubble pa- 
rameter Hq. Together these observations make up dataset 2. 

Throughout, we consider purely Gaussian adiabatic 
scalar perturbations and neglect tensor contributions. We 
assume a universe with a cosmological constant where the 
background cosmology is specified by the following five pa- 
rameters: the physical baryonic matter density Qhh 2 , the 
physical dark matter density Qomh 2 , the ratio of the sound 
horizon to angular diameter distance 6, the optical depth to 
reionisation r and the curvature density f2fc. 

The parameters defining each model are listed in Ta- 
ble [2] together with the ranges of the flat prior imposed on 
them in our Bayesian analysis. 



5 PARAMETER ESTIMATION 
5.1 Power-law parameterisations 

In Fig [2] we plot the 1-D marginal posteriors for the power 
spectrum parameters in each of the power-law models, using 
both dataset 1 and dataset 2. We now on refer to the mean 
of the posterior distribution of each parameter together with 
its 68 % confidence interval. For model 1 (n s ), we see that 
scale-invariant spectrum n s = 1 is ruled out, as expected, 
with high confidence level using either dataset. Indeed, the 
constraints on n s for a flat universe (PL 1) and a curved 
universe (PL 2) using dataset 2 are similar: n s = 0.963 ± 
0.011 and n s = 0.965 ± 0.012 respectively. The inclusion of 
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Figure 2. 1-D marginal posterior distributions for different power 
spectrum parameter models: PL 1 (n s ; top row), PL 2 (n s + Qj.; 
middle row) and PL 3 (n B + n run ; bottom row), using dataset 1 
(solid line) and dataset 2 (dashed line). 



the running parameter (PL 3) displaces and broadens the 
posterior probability, such that n s = 0.944 ± 0.020. 

When the curvature is considered as a free parame- 
ter, dataset 1 selects a closed universe Qk = — 0.033io Q45 
with a Hubble constant Ho = 61. Ol^ 2 , '\ km/s/Mpc. Sim- 
ilarly, for a flat Universe the constraints for the running 
parameter n run = -0.038 ± 0.027 and H = 66.6 ± 3.8 



km/s/Mpc are in good agreement with |Larson et aL (2011 1. 
The inclusion of measurements on different scales (dataset 
2) weakens the geometric degeneracy and tightens signifi- 
cantly the constraints, yielding a mean posterior value for 
the curvature £1^ = —0.0026 ±0.0049 and Hubble parameter 
Ho = 70.4 ± 1.7 km/s/Mpc. On the other hand, we observe 
the best-fit for the running parameter is moved closer to its 
zero value with constraints n run = — 0.018±0.016 

In Fig[3j we plot the reconstructed shape for the power 
spectrum for the simple spectral index in a curved uni- 
verse and the running parameterisation in flat universe, us- 
ing mean values of the posterior distributions found from 
dataset 2. 
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Figure 3. Primordial power spectrum reconstructed using mean 
values of the posterior distributions found from datasct 2, for PL 
model 2 (n s + Q^'t top panel) and PL model 3 (n s + n run ; bottom 
panel), with lcr error bands. 
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Figure 4. Marginalised ID and 2D probability distributions for 
the power spectrum parameters 64, 60 an d A 1 m the LD model, 
using dataset 2. 2D constraints are plotted with lcr and 2cr con- 
fidence contours. 




Figure 5. Marginalised ID probability distributions for the num- 
ber of e-foldings produced during inflation (top) and Hubble con- 
stant (bottom), obtained from the LD model using dataset 2. 



intimately linked to the number of e-folds N during the in- 
flationary epoch through (191. We also find the constraint 
64 [10 -9 ] = — 17.74l|'y| and we note the effect occurring 
through (211: a greater initial curvature (increased I&4J) 



would drive the universe today closer to be flat (decrease 
Ofc — 1|). Finally, the constraint on the scalar field mass is 
M [10- 6 ] = 1.79i° .^ 

The LD model requires a closed universe, then using 
dataset 2 a very small value for the curvature parameter is 
obtained fl k = — 0.43l '2g X 10~ 2 , together with a Hubble 
constant of Ho = 68.9 ± 1.3 km/s/Mpc, both of which are 
compatible with all existing observations. We note, however, 
that many authors have argued on the difficulty to construct 
realistic closed-universe models with the required number of 
e-folds and also that a fine tuning might arise from the initial 
conditions to obtain an homogeneous universe after inflation 



(Ellis etal. 2002 Linde 2003 Uzan et al. 2003 1. In this sense, 



the LD model provides the correct order of magnitude in the 
prediction of the number of e-folds given by TV = 48.1j_4'2 
(see Fig [5b . 

Finally, we plot in Fig. [6] the reconstructed primor- 
dial power spectrum for the LD model together with la 
error bands. It is worth pointing out that the LD spectrum 
automatically contains a clear cut-off at low fc-values and 
an additional fall-off at high wavenumber. A low-fc trunca- 
tion of the primordial spectrum of this type was proposed 
phenomenologically by lEfstathiou ( 2003 1 . This kind of be- 



haviour could be responsible for the small CMB power cur- 
rently observed at low I multipoles. 



5.2 LD model 

We choose to describe the LD spectrum in terms of the pa- 
rameters foo and 64, letting fj, be a derived parameter such 
that the conformal time constraint is fulfilled. The result- 
ing constraints from dataset 2 on the parameters describ- 
ing the LD power spectrum are shown in Fig [4] In par- 
ticular we obtain the constraint bo = 2.45l ;J j7 which is 



6 MODEL SELECTION 

We now investigate which model provides the best descrip- 
tion of the data by performing model selection based on the 
value of the Bayesian evidence. Before performing this pro- 
cess on the real datasets, however, we first test our approach 
by applying it to an idealised dataset. 
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Figure 6. Primordial power spectrum reconstructed from the LD 
model using mean values of the posterior distributions obtained 
from dataset 2, with la error bands. 
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Figure 8. 1-D marginalised parameter constraints for a simulated 
model. The set of parameters used to construct the mock data is 
represented by the vertical lines. The y-axis is relative probability. 
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Table 3. Difference of log-evidences for the different parameter- 
isations relative to simple power-law flat model for the simulated 
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Figure 7. Top (bottom) panel shows a simulated CMB temper- 
ature (polarisation) power spectrum, convolved with chi-square 
noise and limited data up to I of 2000. The model from which it 
was produced is drawn by the blue line. 



6.1 Application to simulated data 

To test our model selection (and parameter estimation) 
method, we implement the LD spectrum in a modified ver- 
sion of the CAMB package (Lewis et al. 2000) to simulate 
the CMB power spectrum predicted by the LD model using 
standard values for the cosmological parameters, i.e. those 
obtained from WMAP7. The predicted temperature and E- 
mode polarisation CMB spectra are plotted as the dashed 
line in Fig. [7| We simulate an idealised process where only 
cosmic variance noise was added to the spectra such that 
each Ci value becomes a random variable drawn from a xiz+i 
distribution with variance 

" '•' (26) 



(AC) 2 = 

The ID marginalised parameter constraints for the LD 
model resulting from the analysis of this simulated dataset 



CMB dataset. 




Model 


A In Z 


PL 1 (n s ) 


0.0 ± 0.4 


PL 2 (n s + n k ) 


-1.0 ± 0.4 


PL 3 (n s + rirun) 


-3.1 ± 0.4 


LD model 


+3.7 ± 0.4 



are shown in Figure [8] We observe that our approach yields 
constraints that appear consistent, within statistical error, 
with the input values used in the simulation, which are in- 
dicated by the vertical lines in each plot. Moreover, these 
results give us an indication of how accurately the param- 
eters of the LD model can be constrained with an optimal 
CMB dataset. In particular, we note that the behaviour of 
the posterior for flk is reflected in that obtained for &4, that 
is because of the link between them through the see-saw 
mechanism. 

The evidence values of the models considered for the 
primordial spectrum are given in Table [3] Based on Jef- 
freys' criterion we observe significant difference in the log- 
evidences for each model. Thus, there is a clear distinction 
between models, with the data clearly indicating a prefer- 
ence, as expected, for the LD model used to generate the 
simulated data. Given the success of model selection on sim- 
ulated data, we now turn to real data with some confidence. 



6.2 Application to real data 

Before performing the Bayesian selection between our dif- 
ferent models for the primordial power spectrum, we will 
momentarily use a frequentist method and compute the 
goodness-of-fit of each of the models using real data. Us- 
ing the best-fit parameters obtained for each model with 
dataset 1, we may calculate Xmim defined as — 21n£ max , at 
the maximum point in each case. These values are given in 
Table 2] from which we see that all three models fit the data 
almost equally as well. 

Turning now to Bayesian model selection, the difference 
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Table 4. Log-likelihood values at the best-fit point for three dif- 
ferent models of the primordial power spectrum using dataset 1. 



Model Wparams -In Cmax 

PL 2 (n s + U k ) 7 7475.3 

PL 3 (n s + rirun) 7 7473.5 

LD model 7 7473.0 



Table 5. Difference of log-evidences for each of our models for 
the primordial power spectrum, relative to the simple power-law 
model. 



Model JVp arams Dataset 1 Dataset 2 

PL 1 (n s ) 6 0.0 ± 0.3 0.0 ± 0.3 

PL 2 (n s + n fe ) 7 -0.5 ± 0.3 -2.6 ± 0.3 

PL 3 (n s + rirun) 7 -0.8 ± 0.3 -1.7 ± 0.3 

LD model 7 +1.2 ± 0.3 -0.9 ± 0.3 



in the log-evidences for each of the models are given in Ta- 
ble [5] as derived from dataset 1 and dataset 2 respectively. 
The log-evidence difference should be interpreted in the con- 
text of the Jeffreys' scale given in Table[T] In the analysis we 
have considered a wide conservative prior for the parameters 
in each model. The parameter estimation results in Section 
[5] showed that the constraints obtained lie well within our 
chosen prior and as such any increase in the prior would not 
include very much more posterior mass within the evidence 
integral. Thus we would not expect reasonable modifications 
of the prior range to alter the results makedly. 

For dataset 1, the similarity in log-evidences and their 
statistical uncertainties make it difficult to obtain definitive 
conclusions. When the evidence values are ranked, however, 
we observe a slight preference in favour of the LD spec- 
trum over the three power-law parameterisations, whereas 
model 3 (n s +n run ) is significantly the most disfavoured. The 
inclusion of additional information on different scales with 
dataset 2 significantly reduces the evidence for the three 7- 
parameter models. We notice that a curved universe with a 
simple power-law, PL 2 (n s + Qk), is strongly disfavoured 
relative to a flat universe with a power-law spectrum (n s ). 
When we perform a comparison of models which contain 
the same number of parameters, we observe the LD model 
is 1.7 units of log-evidence above the PL model 2, which 
is a significant difference, but only 0.8 log-units above PL 
model 3. The ACDM model with a power-law spectrum has 
the largest evidence in this case. We observe from Table [5j 
that from all of the models considered with equal numbers 
of parameters (with either dataset) that the LD model is 
currently preferred. 



can be described in terms of matter-energy content of the 
universe: baryons (f2t>,)i photons (f2 r ), Cold dark matter 
(SIdm) and dark energy (Oa), as well as on its expansion his- 
tory: Hubble time (Ho). The inflationary period is encoded 
in three parameters that have emerged in the construction 
of the appropriate initial conditions. They are related to 
the initial curvature (64), the initial expansion size of the 
universe (bo) and the amplitude of perturbations (fi). Such 
parameters are not predicted a priori, rather we have fit 
them using current observational data. Notice that we have 
the freedom to pick two of them and compute the third pa- 
rameter through the constraint imposed on the conformal 
time. In this paper we have focused on bo and 64 to describe 
the LD spectrum. The comparison of the best-fit primordial 
spectra for our different models is shown on Fig [9] 

We test the LD spectrum by building a simple toy 
model from a chosen CMB spectrum with simulated cosmic 
variance noise. A full Bayesian analysis was performed 
for this toy model and various parameterisations of the 
primordial power spectrum. We have included not only 
the estimation of cosmological and spectral parameters 
but also the Bayesian evidence. For real data, we com- 
pared the x 2 f° r the three different models and concluded 
that the goodness-of-fit for the LD model is about the 
same as the well-studied power-law models. We also 
pointed out that the constraints on the cosmological 
parameters in the LD model are consistent within la error 
bars with the power-law models, as shown in Figs |10| and |ll| 

The Lasenby & Doran model predicts a primordial 
scalar spectrum that at large scales naturally incorporates a 
drop off without the need to parameterise it, while, at small 
scales it automatically incorporates a degree of negative run- 
ning compared with the standard power law parameterisa- 
tion (see e.g. Fig [9] below or Fig. 17 in Lasenby & Doran, 
2005, for a visual comparison). Given that the primordial 
parameters 60 and 64//Z 4 / 3 are well constrained and of order 
unity, no fine tuning is needed for the construction of the 
initial conditions. Moreover, the LD model has the required 
order of e-foldings to obtain realistic models of a closed infla- 
tionary universe. Finally, even though the Bayesian evidence 
from each model showed that it is difficult to make definitive 
conclusions because of the presence of statistical uncertain- 
ties, an important result to emphasise is that for the models 
with the same number of parameters, the preferred model 
to explain current cosmological observations in given by the 
LD model. This may be seen by noting that, for models with 
7 parameters, Table [5] shows that the LD model is preferred 
significantly (by 1.7-2.0 log units of evidence) in all cases 
except relative to n s + n run with dataset 2, for which it is 
preferred by 0.8 log units. 



7 DISCUSSION AND CONCLUSIONS 

The novel closed universe model proposed by Lasenby & Do- 
ran gives rise to a boundary condition on the total available 
conformal time: 

TjM + m = -. 

Here we have split the universe evolution in two main epochs 
that contribute to the amount of conformal time: matter- 
dominated epoch and inflationary period. The matter epoch 
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Figure 9. Primordial and anisotropy power spectrum for differ- 
ent models of the primordial spectrum: n B (dotted line), n e + 
(dash-dot line), n s + n run (dashed line) and LD model (solid line), 
using the best-fit parameter values for dataset 2. 
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Figure 10. Marginalised parameter constraints corresponding to n B (dotted line), n s + £7^ (dash-dot line), n s + n rnn (dashed line) and 
L & D model (solid line), using dataset 1. 
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Figure 11. Marginalised parameter constraints corresponding to n B (dotted line), n s + £7^ (dash-dot line), n s + n rnn (dashed line) and 
L & D model (solid line), using dataset 2. 



